The dielectric response of spherical live cells in suspension: An analytic solution 
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We develop a theoretical framework to describe the dielectric response of live cells in suspensions 
when placed in low external electric fields. The treatment takes into account the presence of the cell's 
membrane and of the charge movement at the membrane's surfaces. For spherical cells suspended in 
aqueous solutions, we give an analytic solution for the dielectric function, which is shown to account 
for the a and /3 plateaus seen in many experimental data. The effect of different physical parameters 
on the dielectric curves is methodically analyzed. 



PACS numbers: 

I. INTRODUCTION 

Dielectric spectroscopy has been successf ully used in 
the past to characterize the biological matter .^^ l ^ l ^ l ^ l ^ l ^ ^^ 
These type of measurements probe the collective dielec- 
tric response of many live cells in suspensions or in tis- 
sues, but by using more or less elaborated models, useful 
information about the state of the individual cell can be 
also extracted.!^ The main advantage of this technique is 
that it is non-invasive, thus the state of the individual 
cells can be monitored without disrupting their natural 
cycle. 

The dielectric response of li ve cel ls is fundamentally 
different from that of dead cells .'^^^ From the dielectric 
point of view, the main difference between the two cases 
is the existence of the membrane potential in live cells. 
The effect of the membrane potential is the accumulation 
of free electric charges at the membrane surfaces. When 
the live cells are placed in time oscillating electric fields, 
these charges move on the surface of the membrane, giv- 
ing rise to extremely high polarizations. Since the mobil- 
ity of these surface charges is relatively small, this effect 
appears only at low frequencies, typically below 10 kHz. 
In this range, the relative dielectric p ermittiv ity of live 
cell suspensions can be as high as IQ^ j^^l^^l IS This phe- 
nomenon is known as the a-relaxation effect. At higher 
frequencies, the a effect disappear and a second interest- 
ing dielectric spectroscopic region appear, which is the /3 
plateau. In this region, the dielectric function of a cell 
suspension is tremendously enhanced by the presence of 
cell's membrane alone. 

In a previous paper,!^ we proposed a model for di- 
electric response of live cells in suspension, which could 
account for both a- and /^-effects. Our focus in that work 
was to give a semi-analytic expression for the dielectric 
function of suspended live cells of arbitrary shape. Be- 
cause of the complexity of such problem, we adopted the 
powerful, but somewhat complicated, spectral method 
introduced by Bergman.l^ In the present work, we work 



out a fully analytic solution of the model proposed in 
Ref.II^, for the case of spherically shaped cells. We hope 
that this analytic solution will provide a new useful tool 
for the scientific community working in this field. 

The independent input parameters for our model are: 
the outer and inner radius of the cell's membrane, ri 
and r2, the dielectric constant and conductivity of the 
medium (eo,cro), of the membrane (ei,cri) and of the in- 
ner cell region (e2, (12), the diffusion constants of surface 
charges accumulated at the outer and inner surfaces of 
the membrane, Di and and the membrane potential 
AV. Excepting 62 and 1^2, all these parameters have 
a very specific effect on the dielectric dispersion curves 
of live cells in suspension, leading us to conclude that 
a fitting of an experimental dispersion curve with the 
present model could provide extremely accurate values 
for all these parameters. 

There is a tremendous amount of theoretical work on 
a-relaxation in colloids. Ref.^^ is one of the earliest work 
that pointed out the fact that the macroscopic dielec- 
tric function of colloidal suspensions is highly dependent, 
in the low frequency range, on the electrical phenomena 
taking place near the surface of colloidal particles. Ref.^^^ 
gave a first semi-quantitative treatment of these effects. 
Refs.^^^^ developed the theoretical foundations for the 
second-layer polarization effect, and the work in the field 
culminated with Ref.^, which gave what was thought 
to be a complete theory of electrophoretic mobility of 
charged colloidal particles. This theory computes the net 
drag force on charged colloidal particles in electrolytes 
and takes into account the deformation of the screen- 
ing cloud via liniarized hydrodynamic equations. How- 
ever, from early comparisons with the experiment, ^^=^ it 
was soon realized that, besides the electrophoretic effect, 
there is another effect, the a-relaxation, which dwarfs 
the first one in most of the cases. The early theory of 
electrophoretic mobilty considered a rigid Stern layer. A 
complete theory will have to relax this assumption. No- 
table efforts in this direction are contained in a string of 
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FIG. 1: (a) Schematic representation of the ion charge ac- 
cumulated around the membrane (blue/red represents posi- 
tive/negative charge distributions), (b) Schematic represen- 
tation of the polarization of this charge when electric fields 
are applied. 
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FIG. 2: In dielectric spectroscopy experiments, a suspension 
of live cells is placed between the plates of a capacitor which 
is subjected to an ac signal, as shown in the left panel. The 
right panel shows our model for live cells. 



papers,^ ^ l ^^ l ^^ l ^^ l which resulted in a fairly complicated 
theory of the polarization of the second layer. A discus- 
sion of latest d evelopments in the field can be found in 

Our model complements these works and we could ar- 
gue that it gives an effective picture of the a-polarization 
of the second layer. The dielectric behavior of live cells 
in suspension at low frequencies and low applied electric 
fields is predominantly determined by the a-relaxation. 
The theoretical model presented here targets specifically 
a-relaxation process, thus allowing us to keep the com- 
plexity and number of model parameters to a mini- 
mum. Other theoretical works specifically addressing 
the dielectric response o f live cells suspensions is con- 
tained in Refs."^^E2l33lMl, Notably, a simple theory of 
a- and /3-effects for spherically symmetric cells, and 
an early attempt for quantifying the /^ -effect for arbi- 
trary geometry. Some later workP^^ give alternative 
approaches to the /^-effect. 



II. MEMBRANE POTENTIAL AND 
ALPHA-RELAXATION 

Live cells contain a large number of negatively charged 
molecules. The inside negative charge attracts positive 
charges from the outside, mainly potassium and sodium 
positive ions. The cells allow most of the potassium 
ions to enter inside, still maintaining an overall negative 
charge, and it keeps most of the sodium ions outside. 
This gives raise to a sharp potential difference across the 
membrane, called the resting membrane potential. Its 
value can be anywhere in a range from 60 mV to a few 
tens over 100 mV. The electric field due to such poten- 
tial differences is enormous. For example, 100 mV over 
a membrane of 10 nm gives an electric field of 10 million 
V/m. 

The charge distribution near the cell's membrane is 
schematically represented in Fig.[l^. As discussed above. 



we have a positive ion distribution outside the membrane 
and a negative ion distribution inside the membrane. The 
exact radial profile of the charge distribution can be com- 
puted by solving the coupled equations describing the 
electrostatics and diffusion of ions near the membrane. 
However, the exact profile is not relevant for the present 
study. What is important is that the distribution is con- 
fined to within a few nano-meters from the surface of the 
membrane. 

When an electric field is applied, the charge distribu- 
tion is deformed as schematically shown in Fig. [TJd. This 
redistribution generates a large dipole moment, which is 
at the origin of the a-relaxation process [cite like a mad 
man]. If the field is time oscillating, the ions will try 
to follow the electric field and dynamically redistribute 
themselves. As we shall see, there is a sharp frequency 
above which the ions can no longer follow the electric 
field. Above this frequency, the polarizability of the cells 
drastically decreases and the a-effect disappear. 

III. THE MODEL 

Our goal is to propose and then solve analytically a 
model for the dielectric response of live cells in suspen- 
sions in the frequency range from 1 Hz to 100 MHz. In 
this range, the electromagnetic field cannot distinguish 
the very fine structure of the cell, but it rather sees an 
effective image of it. It is now generally accepted^^ that 
this effective image is well described by a composite di- 
electric body made of a dielectric shell representing the 
cell membrane [described here by (ei,cri)] and a homo- 
geneous dielectric core [described here by (62,0-2)]. Of 
course, the inside of a cell is very non- homogeneous, but 
this is irrelevant since the field penetrates very little in- 
side the cell. This simple picture of the cell is the start- 
ing point for most of the theoretical work in the field. 
The early work by Maxwell and Wagner, who studied 
the dielectric behavior of spherical dielectric particles in 
suspensions, is probably the most widely used theory to 
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FIG. 3: The interface S separates two electrolytes of different 
dielectric characteristics. A superficial charge distribution p 
is constrained at the interface S. Movement of these charges 
generate the singular current distribution jsing- The diagram 
also shows the usual volume current jVoi in electrolytes and 
the volume Q and its boundary dO. used in the text. 
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FIG. 4: Dispersion curve of eiuj) for: ri = 10 /xm, r2=9.997 /i 
m, (70=0.01, cri=0, (72=0.1, eo=78, ei = 10, £2=80, i:)l = 10"^ 
i:)2=10"^ 71/^1=0.1, 72/i:)2=0.1. The number represent SI 
units. 



describe the dielectric response of cell suspensions.^'' 

On top of this static picture, our model assumes that 
the charge distributions near a cell's membrane can be 
described by effective superficial charge distributions pi 
and p2 at the outer and inner faces of the membrane. 
These superficial distributions of charges are described 
by the following properties: 

1. pi^2 are bound to the faces 111,2 of the membrane, 
so that they cannot leave these surfaces at any time. 

2. The charges are free to move on the faces 2 of 
the membrane. The movement, which is generated by 
gradients in the electric potential and in charge density, 
give rise to singular electric currents: 

Jsing = -7iVs,<^> - AVs.Pi, ^ = 1, 2, (1) 

where 7^ and Di are the electrical conductivities and dif- 
fusion coefficients of the bound charges. 

3. The external electric fields are considered small, so 
that the conductivities 7^ are given with good approxi- 
mation by the the charge distributions poi in the absence 
of any external fields: 

7i = UiPQi, i = 1,2, (2) 

where Ui are the mobility of the bound charges on the 
two membrane's faces. 

From now on, our study focuses on shelled particles 
with the properties described above, which are suspended 
in an electrolyte with dielectric constant eo and conduc- 
tivity ctq. We will use the symbol e* to denote the com- 
plex dielectric defined as e* = e ^ a/juo. We are inter- 
ested in the response of such suspensions when placed 
between the metallic plates of a capacitors like in Fig. |2j 
In the linear regime of small electric fields, the complex 



system will behave like a dielectric material, whose di- 
electric function at given pulsation uo can be computed 
using Lorenz theory: [cite] 

where p is volume fraction occupied by the cells in solu- 
tion and a{(jo) is the frequency dependent polarizability: 

Here, E is the total electric field when a single cell is 
placed in the external, homogeneous, time oscillating 
electric field Eqc'^^^. { ) indicates the average over all 
possible orientations of the cell. 

IV. MATHEMATICAL EQUATIONS OF THE 
MODEL 

We write the equations of the model for arbitrarily 
shaped cells. If we focus on the polarizability 0^(0;), we 
need to consider the problem of a single cell placed in the 
external field EoC^^^. 

The equations governing our model include the Laplace 
equation for the electric potential, V^<^> = 0, the continu- 
ity equations for bound charges, V5]. jsing + dpi/dt = 0, 
together with the boundary conditions at the membrane 
surfaces. The usual boundary conditions at low fre- 
quency have to be modified due to the presence of pi 
at the interfaces. Their new form can be derived from 
the conservation of charges. Indeed, let S be an inter- 
face, with a charge distribution p constrained on it (see 
Fig.[3| and assume that S separates two dielectric media 
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V±. The electrical current flowing near the interface is 
composed by a volume one, given by the usual expression 
= a^E^ and a singular component, jsing flowing on 
the interface. Now, note that, on top of the constrained 
charge distribution p, free charge will accumulate on the 
dielectric interface because of the different conductivities 
cr^. This additional charges behave differently from p, 
because they are free to leave the surface and they don't 
give rise to singular currents at the interface. Denoting 
by r the net superflcial charge distribution, the charge 
conservation for a domain centered on the surface as 
shown in Fig. [3] gives: 



15000 



d 
di 



rdv 



j^o\dS + ^ Jsingdf, (5) 



where T is the contour on the interface given by the in- 
tersection r = dft n S. Using the Maxwell equation, 
r = VD and the divergence theorem, it follows 



on 



{aE^^^D)-dS = - jsingdT. 



(6) 



Noticing that the singular current is the only cause of the 
time variation of the superflcial distribution, variation 



d 
dt 



pdS = (b jsing dT 
nnj: Jr 



(7) 



we arrive at the integral form of the boundary conditions: 

{cTE^^3)-dS=^f pdS. (8) 

In the differential form, the boundary condition takes the 
form: 

nia-E-^l-D-)-nia-E-^ !-&-) = % (9) 

where n represents the normal at the interface (see 
Fig.[3|. 

We can now write the complete set of equations for our 
model: 



V2$ = 0; VsJsing + SpzM = 0. 

E — > Eoex.-p{ju;ot) as |f] oo. 
In the frequency domain: 



V2<l> = 0; Jsing + jc^Pi 
^ Eq as \ f\ oo. 



dt ' 



(10) 



(11) 
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FIG. 5: Dispersion curves of e{oo) for: ri = 10 /xm, r2=9.997 /i 
m, cro=0.001 (green), 0.01 (blue), 0.1 (red), cri=0, ^2=0.01, 
eo=78, £1=10, £2=80, Di^lO'^, D2=10"^ -fi/Di=0.1, 
72/^2=0.1. The number represent SI units. 



V. THE EQUILIBRIUM CONFIGURATION 



In the absence of external electric fields, we can set all 



time derivative in Eq. 10 to zero and obtain: 
V^^ = 0; Vs Jsing = 0. 

<I> ^ as Irl ^ 00. 



0. 



(12) 



The above equations must be complemented by the con- 
dition: 



poi = ±Q 



(13) 



for i = 1,2, respectively. Here, Q is the total charge 
accumulated at the membrane surfaces. 

This system of equations can be easily solved when the 
membrane's conductivity is set to zero, which we will do 
in the rest of the paper. Indeed, the second row of Eq.[l2 
becomes 



aodft^f = and cr2dfi^2 = 0^ 



(14) 



which actually represent trivial Neumann boundary con- 
ditions for the Laplace equation on the outside and inside 
regions of the cell. Consequently, at equilibrium, the po- 
tential is constant inside these regions. Now the equation 
jsing=0 becomes 

Vs, [Ti^tangent - A Vpoi] = - A V ^Oi = 0, (15) 

with a unique solution poi=constant: 

Q 



POi ± 



Area^]. ' 



i = 1,2. 



(16) 
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FIG. 6: Dispersion curves of e(a;) for: ri = 10 /xm, r2=9.997 
^ m, (70 =0.01, (71=0, (72 =0.001 (green), 0.01 (blue), 0.1 
(red), eo=78, ei = 10, £2=80, Dl = 10-^ 02=10'^, j^/Di=0.1, 
72/1^2=0.1. The number represent SI units. 
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FIG. 7: Dispersion curves of e{oo) for: n = 10 /xm, r2=9.99 
/xm (green), 9.995 fim. (blue), 9.997 /xm (red), (7o=0.01, (7i=0, 
(72=0.01, £0=78, £1 = 10, £2=80, Di = 10"^ i:>2 = 10"^ 
71/^1=0.01, 72/1^2=0.01. The number represent SI units. 



Notice that, unlike poi, the total superficial charge at 
the interfaces is not uniformly distributed. Finally, the 
total charge ±(5= J^. TidS is related to the resting mem- 
brane potential by: 



Q = CAV, 



(17) 



where C is the membrane capacitance. 

Since TidS= /^^ PidS, we can state now the main 
conclusion of the section: independent of the shape of 
the cell, the conductivity of the superficial charges pi is 
constant. This constant is controlled by the membrane 



potential: 



7i 



CAV 



Area' 



(18) 



We point out that the diffusion constant D and the con- 
ductivity 7 at the membrane surface are related through 
Einstein relation [/cB=Boltzmann constant]: 



qD = ukuT, 
which, together with Eqs. [2]and[T8) gives: 
/cbT 7i Area^i 



AV = 



q Di C 



(19) 



(20) 



Using the formula for thin capacitors, C=eS/d^ at room 
temperature T = 24^ C, we have: 



AV = 3.00 



7i 



d 



Di ei/ev 



(21) 



where d=ri-r2 is measured in nanometers. The last re- 
lation also shows that the ratios ji/Di and 72/^^2 must 
be the same. 



VI. THE ANALYTIC SOLUTION FOR 
SPHERICAL CELLS 

We will use a single-layer expression for the electro- 
static potential: 



-zEo 

4-K J \f-f'\ ""^^ ' 



(22) 



where we took the oz axis along Eq. Our goal is to solve 
for pi and p2 charge distributions. They are determined 
by the following equations: 



j 7iV£^^ + DiV^^pi = juopi. 
\ e*_iS^^+ - e*S^$7 = Pi. 
On the sphere, we have the following expression: 



d 



d 



1 



Rj sin^ dcP^ 



(23) 



(24) 



where 9 and (j) are the usual spherical coordinates. 

The distributions pi can be expanded in the spherical 
harmonic functions Fim(^,^)- Since the external field is 
homogeneous, only the 1=1 and m=0 term will be actu- 
ally present: 



Aii(r)=Piyio(^), ^2{r^=P2Yio{e). 



(25) 



Using 



47r vi{r, r')Yi^{e, ^1(9', cj>'), (26) 
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FIG. 8: Dispersion curves of e(a;) for: ri=10 /xm, 
r2=9.997 /im, ao^O.Ol, cri=0, (72=0.01, eo=78, ei = 10, 
£2=80, Di^lQ-^, D2=10"^ 7i/i:)i=72/i:>2=0.01 (green), 
0.05 (blue), 0.1 (red). The number represent SI units. 



with 



vi{r, r') 



1 

27Ti^^ 



(27) 



and the ortho-normahzation of the spherical harmonics, 
we obtain: 

^{r)= [-r^o+Pi^i(r,ri)+p2^i(r,r2)]yio(<9). (28) 

(J^o = \J^^o)- Using similar arguments for the pi dis- 
tributions, we write: 

Pi{^ = qiYioiO), P2{^ = q2Yio{e). (29) 

Given that Yio is an eigenfunction of the Laplace operator 
of eigenvalue — 2/i?^, we obtain the following algebraic 
equations for pi , p2 , Qi and q2 : 



(-^1^0 + 


1 r2P2 \ 271 






{-r^Eo + ^ 




+ 1)2^ = 




(el - eS)^o - 


'i+^'or), 1 


2(e*-eS)r2 


= -^1 


(4 - 4)Eo - 


1 ^1 , 




-^2 



(30) 



The solution is: 



Pi 

P2 



— C-B Q^2 ^ 



AC-B 



(31) 



_ 2ri7i + (2Di+ia;r^)(et+2e;) 
2ri7i + (2Di+ja;r2)(e*-e*) 

o ^ 2r27i-2(2Di+ja;r^)(e*-eg)r2/ri 
^ 2ri7i + (2Di+ja;r2)(e*-e*) 



(32) 



c 



r\ 2r272 + (2D2+ja^r^)(e;+2e*) 
r2 2r272 + (2D2+ja;r2)(e*-e*) 



The polarizability of the cells can be computed directly 
from Eq. [sj using E = — and the explicit expression 



of the electrostatic potential, Eq. 28 



3r2 



Pi 
3r2 



P2_ 

3r2 



(33) 



VII. ANALISYS 

In Fig. [4] we report a dielectric dispersion curve gen- 
erated with the present model. The input parameters 
were chosen so that they closely match the experimental 
conditions of live Yeast cells in a buffer solution. The 
volume concentration was chosen p=0.1. In this graph, 
one can clearly distinguishe the a and /3 plateaus, which 
extend from to 100 Hz and 100 to 10^ Hz, respectively. 

In the following, we analyze the effect of variations in 
the different parameters of the model. We start by point- 
ing out that 72 and 62 have very little influence on the 
dielectric properties of the supespension. This confirms 
the assumption made in Ref.^I^ that the charge distribu- 
tion p2 on the inner surface of the membrane plays little 
role in the a-relaxation. This is understandable because 
the electric field penetrates little inside the cell. 

The conductivities of the outer and inner regions of the 
cell, (To and (72, on the other hand, have a very specific 
and similar effects in the P region. From Figs. |5] and [6| 
one can see that <Jo and (72 control the frequency spread of 
the P plateau: the larger ctq or (72 the wider the P plateau. 
For example, by varying (72 from 0.001 to 0.1, which is an 
appropriate range for a inside the cell, we observe a shift 
of the right edge of the P plateau from approximately 10^ 
to 10^. We want to mention that, since the potassium 
ions are taken in by the cells, (72 can be modified by 
changing the potassium concentration of the solution, ctq 
can be modified by, for example, modifying the sodium 
concentration of the solution. 

The thickness of the membrane also has a very specific 
effect on the p plateau. From Fig. [7| one can see that the 
membrane thickness controls the height of the p plateau: 
the smaller the thickness the higher the P plateau. Note 
that the height of the a plateau also changes in Fig. [7[ 
However, the changes in the a plateau are exactly equal 
to the changes in the P plateau, suggesting that the mem- 
brane thickness does not affect the a response of the cells. 

We now focus on Fig. [S] where we fixed d=3 nm and 
let 71/1^1=72/1^2 take the values 0.01, 0.05 and 0.1. 



7 



15000 



10000- 

u 



5000- 








Frequency (Hz) 

FIG. 9: Dispersion curves of e(a;) for: ri = 10 /xm, r2=9.997 
^m, (70 =0.01, (71=0, (72 =0.01, £0=78, ei = 10, £2=80, Di=10~^ 
(green), 10"^ (blue), 10"^ (red), D2 = 10-^ 7i/L>i= 0.1, 
72/^2=0.1. The number represent SI units. 

This imphes the following values for the membrane po- 
tential AV: 9, 45 and 90 mV, respectively. We point 
out that the membrane potential AV of live cells in sus- 
pension can be modified by changing the potassium ion 
concentratioiP^ or by actively blocking or activating the 
ion channels.^ ^ As one can see in Fig. [sj changes in the 
membrane potential have very specific and dramatic ef- 
fects in the alpha region: the larger the membrane po- 
tential, the higher the alpha plateau. For the variations 
in the membrane potential mentioned above, the model 
predicts variations of the alpha plateau of about 8x10^. 

At last, we discuss the effect of the mobility of surface 
charges. Changes in the surface charges mobility also 
have a very specific and dramatic effect in the alpha re- 
gion, as shown in Fig. 9. The larger the mobility, the 
larger the wider the alpha plateau. 

VIII. CONCLUSIONS 

The first goal of our paper was to propose a model 
which can account in a unified way for the dielectric re- 
sponse of live cells in suspensions in both alpha and beta 
regions. The second goal of the paper was to give an 
analytic solution to the model for the simple spherical 
geometry, solution that could be useful for many people 
working in the field. 

Based on this model, we have analyzed the effect of 
different physical parameters on the dielectric dispersion 



curves of live cells in suspension. We found that the con- 
ductivities of the medium and of the intracellular fluid 
control the length of the beta plateau, while the mem- 
brane thickness controls the height of the beta plateau. 
In the alpha region, we found that the membrane po- 
tential controls the height of the alpha plateau while the 
mobility of the surface charges accumulated at the cell's 
membrane controls the length of the alpha plateau. All 
the parameters of the model have distinct influences on 
the dispersion curves, fact that lead us to conclude that 
all the parameters can be accurately obtained by fitting 
experimental data with our model. Thus, the combina- 
tion of experimental dielectroscopy data and our model 
could lead to a methodology for live cell monitoring. 

Our results show that, for a given cell concentra- 
tion and geometry, the low-frequency alpha dielectric 
response correlates with the magnitude of the cellular 
membrane potential. This is potentially very important 
because it enables dielectric spectroscopy to become per- 
haps the only method for monitoring membrane potential 
that is both label free and non- invasive. For example, the 
use of voltage-sensitive dyes is problematic if one wishes 
to simultaneously monitor other physiological processes, 
such as ATP/ADP ratio, that require different fluores- 
cent assays. As a result, if one wishes to study the ef- 
fects of membrane potential on other parameters with 
conventional methods, it is often necessary to use fluo- 
rescent assays on two separate cell populations and then 
to correlate the results with the hope that the two pop- 
ulations are identical. A major advantage of dielectric 
spectroscopy, if further developed, is that one could use 
dielectric response to monitor membrane potential while 
using a fluorescent assay to simulanously monitor another 
parameter on the same cell population. Moreover, dielec- 
tric response appears well suited to flow-cytometry and 
can be readily scaled into multi-electrode (plate reader) 
systems and large scale microfluidic devices. Therefore, 
we believe the method has potential for numerous ap- 
plications, including fundamental research in cell biology 
and biochemistry, pharmaceutical development, and di- 
agnostic methods in medicine. 
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